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ABSTRACT 

The  free  vibrations  of  finite  cylindrical  shells  are  investigated. 

With  the  aid  of  a  number  of  simplifying  assumptions,  a  frequency  equa¬ 
tion  based  on  the  known  characteristic  functions  for  beams  with  any 
combination  of  boundary  conditions  is  obtained.  Experimental  results 
for  frequency  spectra  and  mode  shapes  of  a  cylinder  fixed  on  one  edge 
and  free  on  the  other  are  in  good  agreement  with  both  Rayleigh's  inexten- 
sional  theory  and  the  approximate  frequency  equation.  Structural  damping 
coefficients  obtained  for  the  test  cylinders  are  compared  with  those  of 
previous  investigations. 
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I.  INTRODUCTION 

The  free  vibration  of  a  cylindrical  shell  has  interested  many  investigators. 

In  1894,  Lord  Rayleigh  {Reference  1)  derived  an  approximate  expression  for 
the  natural  frequencies  of  vibration  of  a  cylindrical  shell  based  on  a  separation 
of  the  effects  of  bending  and  stretching.  A  later  treatment  by  Love  (Reference 
2)  resulted  in  a  general  dynamical  theory  of  shells  which  included  both  bending 
and  extensional  deformations.  Love's  equations  were  first  used  by  Flugge 
(Reference  3)  to  obtain  a  cubic  frequency  equation  for  a  simply  supported  cylinder, 
a  result  which  indicated  that  there  were  three  frequencies  for  each  nodal  pattern. 

A  more  detailed  investigation  made  by  Arnold  and  Warburton  (References  4  and  5) 
showed  that  the  three  frequencies  corresponded  to  essentially  radia,  axial  and 
circumferential  vibrations  with  the  radial  vibration  frequency  much  lower  than 
the  other  two.  Their  analysis  also  showed  that  the  natural  frequency  may  de¬ 
crease  as  the  number  of  circumferential  waves  increases,  in  contrast  with  the 
results  of  inextensional  theory.  Arnold  and  Warburton  also  investigated  the 
natural  frequencies  of  cylinders  clamped  at  both  edges,  with  the  use  of  the 
Rayleigh-Ritz  method. 

Recent  investigators  have  concentrated  on  simplifying  the  method  of  analysis 
of  vibrating  cylindrical  shells.  By  means  of  a  number  of  approximations,  Yu 
(Reference  6)  was  able  to  obtain  a  simple  expression  for  the  radial  frequencies 
of  a  clamped  or  simply  supported  cylinder  vibrating  in  a  mode  consisting  of  a 
number  of  circumferential  waves  that  is  large  compared  to  the  number  of  axial 
waves.  Simplified  frequency  equations  were  also  obtained  by  Vlasov  (Reference 
7),  Breslavskii  (Reference  8),  and  Reissner  (Reference  9)  by  neglecting  the  cir¬ 
cumferential  and  axial  inertia  forces  of  the  shell.  Finally,  the  simplifications  of 
Breslavskii  and  Yu  were  combined  by  Rapoport  (Reference  10)  to  yield  frequency 
equations  for  a  shell  with  various  boundary  conditions. 

In  the  present  paper,  a  method  similar  to  Rapoport' s  has  been  used.  An 
experimental  investigation  of  the  frequency  spectra  and  mode  shapes  of  a 
damped-free  cylinder  was  also  performed.  The  experimental  data  are  in  good 
agreement  with  theory.  Structural  damping  was  investigated  as  a  secondary 
part  of  the  experimental  program.  Viscous  damping  coefficients  were  obtained 
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for  each  resonance  point  of  the  supported-free  cylinder  and  tabulated  as  a 
function  of  wave  shape  and  frequency.  The  results  are  compared  with  those 
of  previous  investigations  in  the  Appendix. 

II.  APPROXIMATE  METHOD  OF  ANALYSIS 

The  well-known  Donnell  differential  equation  of  a  circular  cylindrical 
shell  under  an  external  radial  loading  p  can  be  written  as  (Reference  11) 

D  V8 w  +  -2^  -  ?4p  =  0  .  (1) 

R^  3x 


This  equation  can  be  applied  to  vibration  problems  of  cylindrical  shells  when 
we  assume  that  the  circumferential  and  longitudinal  inertia  forces  are  negligi¬ 
ble.  Then  the  external  loading  p  can  be  replaced  by  the  radial  inertia  force 


ph 


On  substituting  this  value  into  Equation  (1)  and  nondimensionalizing  the  resulting 
equation,  we  obtain 


7®w  + 


1 

~~1 


d4  w 

7s* 


+  b 


at1 


-4 
7  w 


=  0 


(2) 


Let  us  assume  w  to  be  of  the  form 


w 


(3) 


Equation  (3)  will  satisfy  Equation  (2)  if  the  coefficients  ^mn  are  the  roots  of 
the  following  equation 


n  = 


(4) 


which  is  an  eighth  order  equation  for  \  as  a  function  of  c,  n  and  0.  Con- 
-  mn 

versely,  we  note  that  if  any  one  of  these  roots  is  known,  then  the  frequency 
parameter  0  is  determined. 
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Since  the  exact  determination  of  X.  is  quite  difficult,  a  number  of 

mn 

simplifying  assumptions  will  be  made  to  obtain  approximate  values.  Let  us 
assume  that  n2  is  large  compared  to  X^.  Then  Equation  (4)  can  be  approxi¬ 
mated  by 


X 


4 

kmn 


c2  (0  n4  -  n8) 


(5) 


For  given  values  of  n,  c,  and  Q  the  right-hand  side  of  Equation  (5)  is  a  constant. 
Therefore,  four  approximate  values  of  kjwnn  are  of  the  form 


mn 


-X 


mn 


iX 


-iX 


mn 


(6) 


The  remaining  four  values  of  X^mn  implied  by  Equation  (3)  are  neglected. 

The  deflection  function  w  can  now  be  written  approximately  as 

w  =  (c.  sin  X  §  +  c,  cos  X  ?  +  c,  sinh  X  § 

'  1  mn  2  mn  3  3  mn 

+  c.  cosh  X  §)  cos  n  8  sin  iDt  (7) 

4  mn  ' 

which  gives  a  longitudinal  deflection  shape  similar  to  that  of  the  vibrating  beam. 
Approximate  values  of  Xmn  are  now  obtained  by  substituting  Equation  (7)  into  the 
appropriate  boundary  condition  equations  for  a  vibrating  beam  and  solving  the 
resulting  determinant.  The  characteristics  roots  Xmn  obtained  by  the  above 
procedure  are  identical  to  the  vibrating  beam  characteristic  roots.  A  tabula¬ 
tion  of  these  values,  for  various  combinations  of  boundary  conditions,  can  be 
found  in  Reference  (12)  and  in  many  text  books  (for  example,  References  13 
and  14).  The  frequency  parameter  Q  for  a  given  n  is  now  obtained  by  substitut¬ 
ing  the  real  values  of  the  beam  characteristic  roots  Xmn  into  Equation  (4)  and 
solving  directly  for  fl. 

Although  the  method  outlined  is  based  on  heuristic  reasoning,  its  justifica¬ 
tion  is  that  the  results  obtained  from  it  are  in  good  agreement  with  experimental 
results.  As  an  initial  check,  results  of  the  approximate  frequency  equation  for  a 
cylinder  with  clamped  ends  are  compared  in  Table  1  with  experimental  frequen¬ 
cies  obtained  by  Koval  (Reference  15)  and  with  the  results  of  another  approximate 
equation  obtained  by  Arnold  and  Warburton  (Reference  5).  In  general,  the  results 


Table  1.  Frequency  Comparison  of  Experimental  Results  with  the  "Approximate 
and  Arnold  and  Warburton  Theories. 
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of  the  present  approximate  theory  are  greater  than  those  obtained  by  Arnold 
and  Warburton  but  are  in  about  as  good  agreement  with  experimental  results. 

From  Table  1  we  see  that  the  agreement  between  the  present  approximate 
theory  and  experiment  becomes  better  as  thr  number  of  circumferential  waves 
increases,  in  accordance  with  the  assumptions  of  the  theory. 

III.  EXPERIMENTAL  INVESTIGATION 

A  series  of  tests  were  performed  on  cylinders  with  one  end  clamped  and 
the  other  free,  as  an  additional  check  of  the  approximate  method  outlined  in  the 
previous  section.  Since  a  detailed  discussion  of  the  test  setup  and  test  pro¬ 
cedure  can  be  found  in  Reference  16,  only  a  brief  outline  will  be  given  here. 

The  two  test  specimens  used  were  made  of  1020  steel  with  dimensions  as 
given  in  Table  1.  The  cylinders  were  formed  over  an  8 -inch  diameter  mandrel 
and  the  seam  was  formed  by  a  butt  weld.  The  cylinders  were  then  spun  to  elimi¬ 
nate  eccentricities.  One  end  of  the  cylinders  was  clamped  in  an  aluminum  plate 
which  contained  a  trough  filled  with  cerrobend,  a  low  melting  point  alloy.  The 
other  end  of  the  cylinder  was  free.  The  cylinder  was  supported  by  a  shaft 
attached  to  the  end  plate  as  shown  in  Figure  1. 

An  electromagnet  was  used  to  excite  the  specimen.  The  test  procedure 
consisted  of  varying  the  frequency  of  the  electromagnet  by  means  of  an  oscillator 
until  a  resonant  frequency  was  reached.  The  frequency  was  accurately  measured 
by  an  electronic  counter.  A  microphone,  which  could  traverse  the  cylinder  axially 
and  circumferentially,  measured  the  response  of  the  shell.  The  microphone  output 
and  the  geometrical  position  of  the  microphone  were  recorded  on  an  x-y  plotter  to 
yield  a  graphical  plot  of  the  longitudinal  and  circumferential  mode  shape  for  a 
given  resonance  frequency.  A  typical  record  is  shown  in  Figure  2. 

IV.  COMPARISON  OF  CLAMPED-FREE  CYLINDER  TEST  RESULTS 

WITH  THEORY 

The  numerical  and  experimental  results  of  the  present  investigation  are 
given  in  Table  2.  For  the  first  longitudinal  mode  of  vibration  (m=l),  the  results 
of  the  approximate  theory  and  of  Rayleigh's  inextensional  theory  are  compared 
with  the  experimental  results  in  Figures  3  and  4.  Numerical  results  of  the 
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Table  2.  Tabulation  of  Experimental  Values  of  Frequency  Spectra, 
Node  Locations,  and  Viscous  Damping  Coefficients. 
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Figure  1.  Over-all  View  of  Test, Setup. 


FREQUENCY  =  3377  CPS 


Page  8 


DISTANCE  AROUND  THE  CIRCUMFERENCE  OF  THE  CYLINDER 
(b)  CIRCUMFERENTIAL  DEFLECTION  SHAPE 
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Figure  4.  Graphical  Comparison  of  Experimental  Results  with  Theory 
for  m  =  1  of  a  Clamped-Free  Shell  (a/h  =  100). 
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approximate  theory  for  both  a  simply  supported-free  cylinder  and  <t  clamped- 
free  cylinder  are  plotted.  All  the  theoretical  curves  are  very  close  to  each 
other  for  the  larger  values  of  n.  It  is  interesting  to  note  that  the  experimental 
values  follow  the  approximate  curve  for  the  clamped-free  cylinder  very  closely 
for  all  values  of  n  for  the  cylinder  with  a  radius  thickness  ratio  of  400.  The 
experimental  results  for  the  cylinder  with  a  radius  thickness  ratio  of  100,  on 
the  other  hand,  fall  between  the  theoretical  curves  for  the  clamped-free  cylinder 
and  the  simply  supported-free  cylinder  for  low  values  of  n  (n  s  4).  The  reason 
for  this  discrepancy  is  suspected  to  be  imperfect  clamping  and  is  being  investi¬ 
gated  further  by  means  of  additional  tests  and  by  a  more  accurate  theory. 

A  comparison,  in  Figures  5  and  6,  of  the  results  of  the  approximate  equa¬ 
tion  with  experimental  results  for  m  a  2  indicates  good  agreement.  The  results 
also  indicate  that  the  value  of  n  at  which  the  minimum  frequency  occurs  depends 
upon  the  axial  wave  length.  As  m  increases,  the  value  of  n  corresponding  to 
the  minimum  frequency  also  increases.  This  phenomenon  was  first  noted  by 
Arnold  and  Warburton  for  cylinders  that  were  simply  supported  or  clamped  at 
both  ends. 

The  position  of  the  experimentally  determined  axial  node  locations  are 
designated  as  in  Table  2.  The  position  of  clamped-free  and  simply  supported 
free  beam  node  points  are  tabulated  in  Table  3.  The  variation  of  the  position  of 
the  experimental  node  points  for  different  values  of  n  is  greatest  at  m-2  and 
decreases  as  the  number  of  axial  waves  increases.  The  tabulated  results  of  the 
average  experimental  axial  node  positions  correspond  fairly  well  with  the  beam 
node  positions  but  do  not  coincide  exactly,  a  result  probably  due  to  the  fact  that 
the  beam  functions  are  not  exact  solutions  of  the  equations  for  the  vibrations  of 
cylindrical  shells. 

V.  CONCLUSIONS 

Experimental  and  theoretical  results  for  clamped-free  and  clamped-clamped 
cylinders  are  in  good  agreement  for  larger  values  of  n  (say,  n  >  4).  It  appears 
therefore,  that  the  approximate  frequency  equation  can  be  used  for  arbitrary 
boundary  conditions  in  this  region.  For  n  >  4  and  m  =  1,  the  Rayleigh  inexten- 
sional  theory  gives  reasonable  results  for  the  clamped-free  cylinders.  The 
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Table  3.  Position  of  Clamped-Free  and  Simply  Supported-Free 
Beam  Node  Points. 


m 

S1 

S2 

S3 

S4 

S5 

i 

0 

- 

- 

- 

- 

2 

0 

0.73/0.79 

- 

- 

- 

3 

0 

0.45/0.51 

0.  85/0.  87 

- 

- 

4 

0 

0. 31/0. 35 

0.  61/0. 65 

0. 89/0. 91 

- 

5 

0 

0. 23/0. 27 

0.47/0.49 

0. 71/0. 73 

0. 93/0. 93 

Key  =  S.  S. -Free/Cl -Free 
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Figure  5.  Graphical  Comparison  of  Experimental  Results  with  Theory 
for  m  >  1  of  a  Clamped- Free  Shell  (a/h  *  400). 


Figure  6.  Graphical  Comparison  of  Experimental  Results  with  Theory 
for  m  >  1  of  a  Clamped-Free  Shell  (a/h  =  100). 
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Figure  6.  Graphical  Comparison  of  Experimental  Results  with  Theory 
for  m  >  1  of  a  Clamped-Free  Shell  (a/h  =  100). 
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experiments  also  indicate  that  although  the  beam  functions  may  not  be  the  true 
deflection  shape,  they  are  near  enough  to  the  true  shape  so  that  when  used  in 
conjunction  with  the  frequency  equation  they  give  a  very  close  approximation 
to  the  experimental  data.  Some  anomalies  in  the  agreement  between  theory 
and  experiment  were  obtained  for  modes  characterized  by  n  <  4.  These  are 
being  investigated  further  at  the  present  time. 
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APPENDIX 

STRUCTURAL  DAMPING 

In  addition  to  the  frequency  spectra,  over-all  structural  damping  coefficients 
were  experimentally  obtained.  These  coefficients  were  obtained  by  passing  the 
response  signal  coming  from  the  microphone  through  an  ac-dc  log  converter  and 
recording  the  output  on  an  oscilloscope  camera.  The  output  was  calibrated  by 
a  db  meter  and  the  decay  curves  obtained  by  disconnecting  the  electromagnet 
from  the  circuit  while  the  cylinder  was  at  resonance.  A  typical  decay  curve 
appears  in  Figure  7. 

The  equation  of  motion. of  the  free  vibration  of  an  elastic  system  (using  the 
same  notation  as  given  in  Reference  17  can  be  written  as: 

a  +  ci  +tii2a  =  0  (c  =  2nf*  g) 
rn  nrn  n  rn  '  n  nB' 

where 

g  structural  damping  coefficient 
f'g  viscous  damping  coefficient 
amplitude  of  nth  mode 
»n  natural  frequency  corresponding  to 

Tabulated  values  of  f'g  as  a  function  of  mode  shape  and  frequency  are  given  in 
Table  2.  A  plot  of  the  "viscous"  damping  coefficient  as  a  function  of  the  number 
of  circumferential  waves  is  shown  in  Figure  8.  The  results  indicate  a  large 
scatter  with  f'g  s  2  being  an  average  value  for  both  cylinders.  It  is  interesting 
to  note  that  an  average  value  of  two  for  the  viscous  damping  coefficient  is  clo^e 
to  that  found  by  the  author  in  Reference  16  for  a  steel  cylinder  clamped  on  both 
ends.  Fung,  Sechler  and  Kaplan  obtained  an  average  value  of  f'g  =  6  for  a  set 
of  aluminum  cylinders.  Figure  8  also  shows  that  for  n  <  8,  the  m  =  1  and  m  =  2 
modes  have  an  average  value  of  one  for  f'g  and  the  m  =  4  and  m  =  5  modes  an 
average  value  of  three.  The  m  =  3  modes  seem  to  fluctuate  between  one  and 
three.  The  peaking  effect  at  a  unique  frequency  found  by  the  author  in  Reference 
16  did  not  occur  in  the  present  Investigation.  It  is  clear  that  much  work  remains 
to  be  done  before  the  structural  damping  phenomenon  is  completely  understood. 
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